Debye mass and heavy quark potential in a PNJL quark plasma 
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We calculate the Debye mass for the screening of the heavy quark potential in a 
plasma of massless quarks coupled to the temporal gluon background governed by 
the Polyakov loop potential within the PNJL model in RPA approximation. We give 
a physical motivation for a recent phenomenological fit of lattice data by applying 
the calculated Debye mass with its suppression in the confined phase due to the 



Polyakov-loop to a description of the temperature dependence of the singlet free 



■ energy for QCD with a heavy quark pair at infinite separation. We compare the 

O • result to lattice data. 
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I. INTRODUCTION 

The proposal of Matsui and Satz [1] that color screening of heavy quarks in a deconfined 
medium should lead to a dissociation of the J/i/j bound state and whence could be a clear 
signal of quark-gluon plasma (QGP) formation in heavy-ion collision (HIC) experiments 
arose to a wide field of research both in experimental as well as theoretical directions j^J. 
Following the original idea the simplest picture is to consider a static quark-antiquark (QQ) 
probe in a color plasma where all medium effects are included in two body interaction poten- 
tial Vr{r). Then solving the QQ T-matrix |3| for such a system (or equivalently Schrodinger 
equation [4]) one can obtain a temperature modification of the hadron spectra. For obvious 
reasons this description is unrealistic (e.g., in HIC experiments the medium always has col- 
lective flow) but nevertheless it can provide physically reliable insights. But what potential 
should we use? One possibility of approaching this problem is to use nonperturbative field 
theoretical methods to motivate an ansatz for the heavy quark potential and then fit its 
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parameters to the lattice data. For example in ref. {3, 5| it has been shown that including 
effects of a dimension-two gluon condensate one gets quarkonium spectral functions in good 
agreement with lattice data. 

In this note we go similar way by investigating Debye screening effects in modifying the 
vacuum Cornell potential (one-gluon exchange plus linear term) when going to a finite tem- 
perature medium. The screening medium is composed of massless quarks coupled to the 
homogeneous, temporal gluon background which is a very crude approximation of confine- 
ment. This essentially constitutes a version of the PNJL model [8|, |9j which was found to 
reproduce lattice thermodynamics although it is theoretically still rather unsatisfactory. In 
this contribution we discuss the similarities of the temperature dependence of the Debye 
mass within such a model and that obtained from a fit of lattice QCD heavy quark free 
energy with the effective interaction of Ref. [3] . This observation holds promise for a better 
microscopic understanding of the physical effects determining the behaviour of the heavy 
quark- ant iquark potential at finite temperature. 



II. SCREENING IN THE GLUON BACKGROUND 
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In Ref. p] a model was considered where the vacuum potential was screened by quark 
matter polarization decribed by quark-antiquark loop integrals where the internal lines were 
coupled to a temporal gluon background field. Specifically, when the static interaction po- 
tential is given as V(q), q 2 = |q| 2 , the statically screened potential is given by a resummation 
of one-particle irreducible diagrams ("bubble" resummation = RPA) 

Kc(g) = V(g)/[l-n 00 (0;q)/g 2 ] , (1) 

where the longitudinal polarization function in the finite T case can be calculated within 
thermal field theory as 

n o(^; q) = 9 2 T [ 7^~Ti Trh S$(iu n ; p)j S^(iu n - iu t ; p - q)] . (2) 

Here 00 1 = 2nlT are the bosonic and u n = (2n+l)irT are the fermionic Matsubara frequencies 
of the imaginary-time formalism. The symbol Tr stands for traces in color, flavor and Dirac 
spaces. S$ is the propagator of a massless fermion coupled to the homogeneous static gluon 
background field (fs. Its inverse is given by js, 9] 

^(p; w„) = 7 ■ P + y iu n ~ 70-W3 , (3) 



is, 



where <p 3 is related to the Polyakov loop variable defined by 

$(T) = iTr c (e^) = ^(1 + 2cos(/3y? 3 )) ■ 

The physics of $(T) is governed by the temperature-dependent Polyakov loop potential 
W($), which is fitted to describe the lattice data for the pressure of the pure glue system js}. 
After performing the color-, flavor- and Dirac traces and making the fermionic Matsubara 
summation, we obtain in the static, long wavelength limit 

Iloo(q) = r dpp^ = -2 g 2 T 2 I(<S>) = -ml(T) , (4) 

JO Op 

where m£>(T) is the Debye mass, the number of degrees of freedom is N c = 3, Nf = 2 and 
f<s>(p) is the quark distribution function {{J. In comparison to the free fermion case 10, ll| 
the coupling to the Polyakov loop variable $(T) gives rise to a modification of the Debye 
mass, given by the integral 

/($) = —/ dxx — — — . (5) 

y J 7T 2 Jo 1 + 3$(1 + e- x )zz x + e~ 3x v ; 

The temperature dependence of $(T) is taken from Ref. 12|. In the limit of deconfinement 

($ = 1), the case of a massless quark gas is obtained (1(1) = 1), while for confinement 

($ = 0) one finds that 1(0) = 1/9. The temperature dependence of the resulting Debye mass 

is shown in Fig. [1] and as expected from the very beginning is much lower, comparing to the 

free, massless case, in the confined and transition region (with the critical temperature T c w 

200 MeV). For temperatures T » T c the free gas behavior is reproduced. Adopting the field 

theoretical approach of Ref. that the color singlet free energy at finite temperature 

is driven by the screened, nonperturbative gluon propagator and that Debye masses in 

Coulomb and stringy sectors are different, we get (a s = g 2 /An) 

, . 4 ( e - mDr ml m 2 r ~ r \ 

F Qo (r) = — a s \ ^3- + -J^e" mDr + m D , 6 

QQV ; 3 V r 2m D 2m D J ' V ; 

where the constant (r- independent) term is a homogeneous mean field contribution and 

takes into account the one-particle self-energy effects We take m 2 G = 0.631 GeV 2 

which gives string tension a = 2m 2 G a s j?> = 0.198 GeV 2 . We notice that our calculated 

Debye mass resembles qualitatively the behavior resulting from a recent fit for the stringy 

term in Ref. (3] with the difference that our result is smaller than that of the fit by abut a 

factor two. Therefore we proceed by plugging it in for m~D in equation ((6]), and if we define 

Foo(T) = Fqq(t -> oo) = \a s - m D ) . (7) 
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Then we see that F^T — > oo) ~ — T, and that it diverges for small temperatures - faster 
than in the free fermion case due to gluon suppression of tHd - which is expected in quenched 
approximation but not in QCD with dynamical quarks which would make Foo finite due to so 
called "string breaking mechanism" . Fitting lattice data for Foo(T) 2j for high temperatures 
results in a coulombic Debye mass mo = 0.12 T which is much lower than the perturbative 
one which is about 3.143 T for two flavors and a s = 0.471. We plot the resulting free energies 
in Fig. [TJ What concerns F^lT) we see that our results meet quite well the lattice data in a 
quite a wide temperature range starting from the transition region up to high temperatures; 
the generic shape is reproduced in the transition region, while we see a discrepancy for low 
temperatures. Our result is smaller than the lattice (especially for low T) which is a remnant 
of our lower Debye mass than the one of the fit by Riek and Rapp 2J- What concerns the 
free energy for small distances, it is temperature independent and matches T = potential, 
while r rj 0.25 fm marks the onset of temperature effects. The temperature dependence of 
binding energies obtained with this approach suggests f| that charmonium (bottomonium) 
dissociation takes place at T « 200 MeV (T w 250 MeV). 



III. CONCLUSIONS 



It is clearly shown that even a crude inclusion of residual effects from the gluon sector 
results in a remarkable improvement. One of the obvious extensions of this calculation is to 
take into account chiral symmetry effects on the internal quark lines, governed in this model 
by the local NJL interaction [9]. In this way one could in principle study the effects of the 
chiral transition on J/ip dissociation. 
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♦ Lattice data [hep-lat/0506019] 




1.2 

> 

a) 0.9 
O 

"o0.6 

s 

0.3 



with Polyakov loop 
without Polyakov loop 




FIG. 1: Lower panel: Temperature dependence of the Debye mass with (solid line) and without 
(dashed line) inclusion of Polyakov loop. Calculated for a s = 0.471. Upper panel: Temperature 
dependence of F\(r = oo,T) with (solid line) and without (dashed line) inclusion of Polyakov loop, 



compared to lattice data from Ref. 
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